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ABSTRACT 



> 

<N 

. We have acquired between 127 and 150 h of time-resolved multicolour photometry 

f*^ ' for each of the three j3 Cephei stars IL Vel, V433 Car and KZ Mus over a time span of 

£T} , four months from two observatories. All three objects are multiperiodic with at least 

f^*) ■ three modes of pulsation. 

^5 ' Mode identification from the relative colour amplitudes is performed. We obtain 

unambiguous results for the two highest-amplitude modes of IL Vel (both are t = 1) 
and the three strongest modes of KZ Mus (I = 2,0 and 1), but none for V433 Car. 
Spectroscopy shows the latter star to be a fast rotator (vs'mi = 240km/s), whereas 
the other two have moderate wsini (65 and 47 km/s, respectively). 

We performed model calculations with the Warsaw-New Jersey stellar evolution 
and pulsation code. We find that IL Vel is an object of about 12 Mq in the second 
half of its main sequence evolutionary track. Its two dipole modes are most likely 
rotationally split components of the mode originating as p\ on the ZAMS; one of 
rS \ these modes is m = 0. V433 Car is suggested to be an unevolved 13 Mq star just 

■ entering the (3 Cephei instability strip. KZ Mus seems less massive (w 12.7Mq) and 

somewhat more evolved, and its radial mode is probably the fundamental one. In this 
case its quadrupole mode would be the one originating as g±, and its dipole mode 
would be pi . 

Two of our photometric comparison stars also turned out to be variable. HD 90434 
is probably a new slowly pulsating B star whose dominant mode is a dipole, whereas 
the variability of HD 89768 seems to be due to binarity. 

It is suggested that mode identification of slowly rotating (3 Cephei stars based 
on photometric colour amplitudes is reliable; we estimate that a relative accuracy of 
3% in the amplitudes is sufficient for unambiguous identifications. Due to the good 
agreement of our theoretical and observational results we conclude that the prospects 
for asteroseismology of multiperiodic slowly rotating [3 Cephei star are good. 

Key words: stars: variables: other - stars: early-type - stars: oscillations - stars: in- 
dividual: IL Vel - stars: individual: V433 Car - stars: individual: KZ Mus - techniques: 
photometric 



1 INTRODUCTION 

The j3 Cephei stars are a group of early B-type stars of lu- 
minosity classes III-V that pulsate in pressure (p) and grav- 
ity (g) modes of low radial overtone. Their pulsational be- 
haviour is thus similar to that of the A/F-type 5 Scuti stars. 



Both classes of variable seem suitable for asteroseismologi- 
cal investigations: their pulsational frequencies may be used 
to sound their interiors as the associated modes penetrate 
deeply into the star. 

As exciting as this possibility is in theory, practice has 
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shown that several obstacles still need to be overcome before 
precision asteroseismology of /3 Cephei or 8 Scuti stars can 
be done. On the theoretical side an adequate description 
of higher-order rotational effects is being worked on (e.g. 
Soufi, Goupil & Dziembowski 1998), and mode coupling is 
expected to have major effects on the pulsational eigenspec- 
tra of fast rotators (see Daszyriska-Daszkiewicz et al. 2002). 
Observationally, the main problem is the detection and iden- 
tification of as many pulsation modes as possible because a 
complete set of mode spectra has not even nearly been de- 
tected for any of these pulsators. 

It may be suspected that the S Scuti stars offer better 
possibilities for asteroseismology than ft Cephei stars do, as 
they generally have many more pulsation modes excited to 
observable amplitudes. This result is in agreement with the- 
oretical investigations of mode excitation (e.g. Pamyatnykh 
2003). However, the problem lies with mode identification, 
that is, the match of the observed frequencies with the cor- 
responding pulsation mode defined by the three parameters 
k, the radial overtone of the mode, the spherical degree I 
and the azimuthal order m. 

To accomplish such a match in the presence of in- 
complete observed pulsation spectra, it is necessary to use 
mode identification methods. However, the surface convec- 
tion zones of 8 Scuti stars are expected to affect some of 
these methods, like the commonly used photometric method 
that utilises and amplitude ratios and/or phase shifts of the 
light curves in different filters, adversely (see Balona & Ev- 
ers 1999), and it is not known if they can be trusted. The 
(3 Cephei stars do not have surface convection zones and 
therefore the simple photometric method may be applicable 
with more confidence. 

We have therefore chosen three /3 Cephei stars whose 
literature data qualify them as interesting candidates for 
in-depth studies aiming at the detection of many pulsation 
frequencies and at their mode identification by photometric 
means. All these objects have been reported as multiperiodic 
variables before and photometric mode identifications have 
been attempted. However, in all cases the amount of data 
available was small enough to give us hope for a significant 
improvement over previous results if we can obtain ^ 100 
h of observation. 

The most comprehensive photometric study of IL Vel 
(HD 80383) was performed by Heynderickx & Haug (1994). 
They found four frequencies in their UBV light curves, three 
of which formed an equally spaced triplet, suggestive of ro- 
tational m-mode splitting. This frequency solution did how- 
ever not fit additional data in the Walraven system within 
the precision of the measurements. Heynderickx, Waelkens 
& Smeyers (1994) used the Walraven data for mode identifi- 
cation and suggested that all modes are of degree £ = and 
t = \. 

The variability of V433 Car (HD 90288) was discov- 
ered by Lampens (1988) and followed up by Heynderickx 
(1992). The frequency analysis of all the Geneva and Wal- 
raven data available to Heynderickx (1992) allowed him to 
disentangle four frequencies in the light variations of V433 
Car. The mode identification by Heynderickx et al. (1994) 
was somewhat ambiguous, but suggested a mixture of (pos- 
sibly) radial and nonradial modes with spherical degrees i 
up to 4. 

Our third target star, KZ Mus (HD 109885) was a dis- 



covery of the HIPPARCOS satellite (ESA 1997) reported by 
Waelkens et al. (1998). Aerts (2000) obtained new Geneva 
photometry for the star and performed both a frequency 
analysis of all data and mode identification. She suggested 
that the star pulsates with at least three frequencies, where 
the second strongest one appeared to be radial. 



2 OBSERVATIONS 
2.1 Photometry 

We have studied IL Vel, V433 Car and KZ Mus from both 
the Siding Spring (SSO, Australia) and South African Astro- 
nomical (SAAO) Observatories. We used photoelectric pho- 
tometers attached to the 0.6-m telescope at SSO and the 
0.5-m and 0.75-m telescopes at SAAO. 

The most valuable part of the energy distribution of a 
(3 Cephei star for mode identification by means of photom- 
etry is in the blue (A < 4200 A). Therefore, the Geneva and 
(particularly) Walraven systems are well suited for such a 
purpose as they have several filters in this wavelength do- 
main. Unfortunately, these are not widely available. Owing 
to the faintness of our variables (V = 8.1 — 9.1) we decided 
to utilise the standard Johnson UBV system, supplemented 
by Str0mgren v whenever possible. 

The periods of the known pulsations of our three targets 
are all longer than 150 minutes, and the stars are all located 
in an area covering some 30° on the sky. Consequently, we 
chose one local comparison star for each variable (the B5 
IV star HD 79670 for IL Vel, the B9 IV/V star HD 90434 
that was later replaced by the Al IV/V star HD 89768 for 
V433 Car and the B9 IV star HD 109082 for KZ Mus; all 
spectral types are from Houk & Cowley 1975) , and our ob- 
serving sequence included all six stars whenever reachable. 
The resulting cycle time of about 25 minutes (or less) per 
measurement of each variable sampled the light curves of all 
stars adequately. 

Another aspect of the pulsations of (3 Cephei stars re- 
quires consideration when planning an observational effort 
such as ours. In many cases, closely spaced pulsation fre- 
quencies are present. As mentioned above Heynderickx & 
Haug (1994) reported just that for our target IL Vel. Con- 
sequently, a sufficient time baseline must be spanned by the 
observations (at least 2 months for IL Vel) for such frequen- 
cies to be resolved and the measurements should be dis- 
tributed evenly over the observation period to avoid aliasing 
problems. 

We have therefore also observed these objects during a 
multisite campaign for the 8 Scuti star FG Vir (Breger et 
al., in preparation) when this star was not reachable. This 
only resulted in short runs up to two hours that however 
covered important gaps in the total coverage which made 
them quite valuable. The main body of our final data set 
spanned 94 days in February to May 2002, during which data 
was obtained on 52 nights. A total of 146 hr of measurement 
of IL Vel, 150 hr of V433 Car and 127 hr of KZ Mus was 
acquired. 

Data reduction was performed in the standard way for 
differential photoelectric time-series photometry. First, we 
corrected the data for coincidence losses, followed by sky 
subtraction and by extinction determination using the time 
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Figure 1. V- filter light curves from all our measurements of IL Vel (filled circles), V433 Car (open circles) and KZ Mus (asterisks). Note 
the shorter periods and much lower amplitude of V433 Car. Multipcriodicity of all three stars is evident. 



series of the constant stars in our ensemble. After the extinc- 
tion correction, the differential light curves of the variables 
with respect to their local comparison stars and that of the 
three comparison with respect to the others were calculated 
and the times of measurement were converted to Heliocentric 
Julian Date (HJD). The light curves of the three /3 Cephei 
stars in the V filter are shown in Fig. 1. 

As part of the reductions were already performed during 
the observational campaign, it was noticed that the com- 
parison star HD 90434 was variable. Consequently, it was 
replaced by HD 89768. As it turned out during the final 



reductions after completion of the observations, this star is 
also variable. We will analyse these objects later, but for 
the time being let it suffice to say that we performed a fre- 
quency analysis of their differential light curves with respect 
to the constant comparison stars for the other targets. We 
then subtracted synthetic light curves built from the results 
of this frequency search from the differential measurements 
of V433 Car so the variability of its comparison stars will 
not affect the frequency analysis of this star. This procedure 
is valid as the time scales of the variability of HD 90434 and 
HD 89768 is much longer than that of V433 Car. 
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Table 1. Results of our determinations of the radial and pro- 
jected rotational velocities. 



Star 


vsini [km/s] 


< Vr > [km/s] 


IL Vel 


65 ± 3 


+ 19.0 


V433 Car 


240 ± 10 


+4.0 


KZ Mus 


47 ± 3 


-61.6 



2.2 Spectroscopy 

With the main aim of determining the projected rotational 
velocities of IL Vel, V433 Car and KZ Mus, high-resolution 
spectra of the three f3 Cephei stars were obtained on the 
night of 4/5 May 2002 with the 1.9-m telescope at SAAO. 
We used the GIRAFFE echelle fibre-fed spectrograph attached 
to the Cassegrain focus of this telescope. The GIRAFFE spec- 
trograph has a resolving power of about 32 000, giving a 
resolution of 0.06 - 0.09 A per pixel. Exposure times were 
1500 seconds for IL Vel and KZ Mus as well as 600 seconds 
for V433 Car for 7S/N ratios between 27 and 35. A Th-Ar 
arc lamp was used for wavelength calibration. 

The spectra were normalised to the continuum by using 
an unbroadened synthetic spectrum with T e fr = 23000 K, 
logp = 4.00 as a template, using the SPECTRUM code (Gray 
& Corbally 1994). A running median of each echelle order 
was divided by the corresponding section of the synthetic 
spectrum and taken to represent the response function of 
the instrument. A polynomial of degree 5 was fitted to the 
response function and used to correct the observed spec- 
trum. 

For each order, the observed rectified spectrum was cor- 
related with the corresponding section of the synthetic spec- 
trum after removing the unit continuum, effectively result- 
ing in the mean line profile with the continuum removed. A 
quadratic was then fitted to the correlation function, and its 
maximum adopted as the "radial velocity" for each order. 
The mean radial velocity from all the orders is also obtained. 

The correlation profiles can be co-added to form what 
is essentially a mean line profile. The projected rotational 
velocity, v sini, can be determined by fitting a model profile 
of a rotating star. The projected rotational velocity is ad- 
justed until a best fit is obtained to the observed profile. We 
show these fits in Fig. 2; the results from our spectroscopic 
analysis are summarised in Table 1. 



3 FREQUENCY ANALYSIS OF THE 
PHOTOMETRIC TIME SERIES 

Our frequency analysis was performed with the pro- 
gram PERIOD 98 (Sperl 1998). This package applies single- 
frequency power spectrum analysis and simultaneous multi- 
frequency sine-wave fitting. One of its many advanced fea- 
tures is its capability to fix dependent signal frequencies 
to certain values, e.g. to sum and difference terms, and to 
perform simultaneous nonlinear least squares fits with such 
fixed frequencies. 

Our strategy for the frequency analysis starts with the 
calculation of the spectral window of the data. It is com- 
puted as the Fourier amplitude spectrum of a single noise- 
free sine wave with the frequency of highest amplitude in 




Figure 2. The average Helium line profile for our three vari- 
ables (full line) overlaid by a model atmosphere fit (dashed lines) 
to determine the projected rotational velocities. The line profile 
asymmetries are likely caused by nonradial stellar pulsation. 



our data. In this way, the reflection of the spectral window 
pattern at zero frequency, which may not be negligible as 
our observing nights often did not sample a full cycle of the 
light variations, can be evaluated. 

We continue by computing amplitude spectra of our 
data as well as those of residual light curves after the pre- 
viously identified periodicities had been removed using the 
multi-periodic fitting algorithm of PERIOD 98. We continue 
this process until no significant peaks are left in the resid- 
ual amplitude spectrum. We consider an independent peak 
significant if it exceeds four times the local noise level in the 
amplitude spectrum, following Breger et al. (1993). Combi- 
nation signals only require S/N > 3.5 to be regarded as 
significant. Experience has shown that this criterion is both 
reliable and conservative. 

Similar analyses were performed for all the four filters 
used. We note that, within the uncertainties of frequency 
determination, these detected frequencies were the same for 
all stars in every filter. This confirms that alias ambiguities 
do not affect our analysis. As our final values for the fre- 
quencies we adopted the mean values from the data in all 
four filters weighted by the associated signal-to-noise ratios. 
The deviations of the optimum frequencies in the individual 
filters from this mean is taken as our error estimate on the 
frequency determination. 



3.1 KZ Mus 

As the frequency analysis of this star was straightforward, 
it is the first that we chose to describe. The spectral win- 
dow, amplitude spectra, and prewhitened versions thereof 
are shown in Fig. 3. We use the data obtained in the B filter 
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Figure 3. Spectral window and amplitude spectra of our B filter 
data of KZ Mus. Some prewhitening steps of the detected fre- 
quencies arc shown in consecutive panels. The lowest panel shows 
the amplitude spectrum of the combined UvBV residuals after all 
significant signals have been removed from the measurements. 



for display as they have best S/N ratio. Due to our two-site 
coverage the spectral window is sufficiently clean so that 
alias ambiguities do not represent a problem for the analy- 
sis. 

The first three signals in the periodogram are easy to 
detect. After prewhitening these, a more complicated resid- 
ual amplitude spectrum (fifth panel of Fig. 3), with several 



peaks standing out, remains. A closer look reveals that two 
of them are the combination sum and difference frequencies 
of the two strongest modes, whereas the remaining tall peak 
corresponds to an independent signal. Prewhitening all the 
six variations detected so far, the 2f-harmonic of the second 
strongest mode also becomes notable. 

After this signal is included in our multifrequency light- 
curve fit, no significant peaks remain in the residual ampli- 
tude spectrum of the B filter data alone, and neither are 
there any present in the measurements in any other filter. 
However, we can attempt to increase the signal to noise in 
our analysis by adding the data in all filters. Obviously, this 
must not be done for the original light curves, as differ- 
ent pulsation modes have different relative amplitudes in 
the different wavelength passbands (that are a diagnostic 
for mode identification!). It is however safe to combine the 
residual light curves as we are working at low S/N ratio, 
i.e. the intrinsic differences in mode amplitude between the 
wavebands are small compared to the noise in the data. 

To combine the different residual data sets we have 
therefore subtracted a multifrequency solution from all the 
individual filters' observations, using exactly the same fre- 
quencies (determined as described in the previous section) 
for each waveband. We then roughly scaled the residuals to 
B amplitude by multiplying the U filter residuals by 0.8 and 
the V filter residuals by 1.1. We did not scale the v filter 
data as the signals in it were found to have essentially the 
same amplitude as in B. The scale factors are approximate 
average amplitude ratios of the modes previously detected. 

All these residual light curves were then added together, 
averaged, and the periodogram was computed out to the 
Nyquist frequency (lowest panel of Fig. 3); we find no more 
significant signals. Interestingly, we do note two residual 
mounds of amplitude, one in the frequency region where 
the pulsational signals were found and another one where 
combination frequencies would be expected. In fact, the two 
highest peaks in the latter region are located at exact sums 
of two known frequencies but they are not formally signifi- 
cant. 

We can therefore claim the detection of seven signals in 
the light curves of KZ Mus; our final frequency solution is 
listed in Table 2. The rms scatter per single data point in 
the combined and averaged UvBV residuals is 3.0 mmag. 

3.1.1 Re-analysis of published data 

Our frequency solution in Table 2 differs somewhat from the 
results of Aerts (2000), who needed to base her frequency 
analysis on the HIPPARCOS data of the star. The first two 
signals detected by her agree with ours, but the third one is 
different; the remaining signals are below her detection level. 
Two possible reasons for the disagreement can be imagined, 
the more exciting one being intrinsic amplitude variability 
of the star and the other one being an alias problem in the 
poorly sampled HIPPARCOS photometry. 

Consequently, we have re-analysed the time series of 
KZ Mus acquired by the HIPPARCOS satellite. The resid- 
ual amplitude spectrum after prewhitening the two consis- 
tent main frequencies contains the third signal claimed by 
Aerts (2000) as the tallest peak. However, our third signal 
is also visible and of similar strength. As our data are more 
extensive and of better quality than the HIPPARCOS pho- 
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Table 2. Multifrequency solution for our KZ Mus data. Frequencies arc mean values from results of all four filters weighted by the S/N 
ratio of the signal. Errors on the amplitudes were calculated with the formulae of Montgomery & O'Donoghuc (1999). The S/N ratio 
quoted is the average of the values in the individual filters. 



ID 


Frequency 




Amplitude 




mean S/N 




(c/d) 


U (mmag) 


B (mmag) 


V (mmag) 


v (mmag) 




h 


5.86384 ± 0.00006 


45.6 ± 0.3 


40.9 ± 0.3 


38.5 ± 0.3 


41.2 ± 0.4 


89.4 


h 


5.95026 ± 0.00011 


28.8 ± 0.3 


19.6 ± 0.3 


16.4 ± 0.3 


20.6 ± 0.4 


45.6 


h 


6.1874 ± 0.0005 


14.3 ± 0.3 


11.0 ± 0.3 


10.5 ± 0.3 


11.8 ± 0.4 


25.7 


h 


5.7090 ± 0.0024 


3.4 ± 0.3 


3.2 ± 0.3 


2.4 ± 0.3 


3.1 ± 0.4 


6.6 


h-h 


0.08642 


4.3 ± 0.3 


2.8 ± 0.3 


3.1 ± 0.3 


3.1 ± 0.4 


6.3 


h + h 


11.81410 


2.5 ± 0.3 


2.2 ± 0.3 


2.5 ± 0.3 


2.6 ± 0.4 


5.7 


2/2 


11.90052 


2.6 ± 0.3 


1.9 ± 0.3 


2.0 ± 0.3 


2.6 ± 0.4 


5.2 



tometry, we suggest that the third signal found by Aerts 
(2000) is an unfortunate artifact of the poor sampling of the 
HIPPARCOS measurements. 

3.2 V433 Car 

The upper panel of Fig. 4 shows that the spectral window 
function for our measurements of this star is also quite good. 
The amplitude spectrum is dominated by two peaks of sim- 
ilar strength that seem to imply an aliasing problem, but 
thanks to the long time base of our data they are well re- 
solved and do not interfere with each other; they are simply 
two independent signals of similar amplitude. 

Prewhitening of these variations results in the detection 
of a third frequency and the next step reveals another signif- 
icant peak. After these four periodicities have been removed, 
no obvious new candidate frequency can be discerned. We 
therefore combined the residual data as we did for KZ Mus, 
but used somewhat different scale factors that seemed more 
appropriate for this star. We multiplied the U residuals by 
0.83 and the ones in V by 1.05 to obtain approximate B 
amplitudes. 

The amplitude spectrum of the residuals thus combined 
is shown in the lowest panel of Fig. 4, again out to the 
Nyquist frequency. The highest peak corresponds to a com- 
bination frequency difference and as it exceeds the noise level 
more than 3.5 times, we adopt it as an additional term for 
our frequency solution. We would also like to point out the 
peak near 7 c/d in the lowest panel of Fig. 4. It is in fact 
a close doublet at 6.924 and 6.968 c/d that is present in 
the residuals of all individual filters, but as it still does not 
exceed S/N > 4, we cannot be sure about its reality. We 
therefore adopt a five- frequency solution for our light curves 
that we list in Table 3. 

The rms scatter of the combined residual light curves 
of V433 Car is 3.9 mmag per point. This is 30% higher than 
for KZ Mus. Likely reasons for the increased noise are un- 
detected further intrinsic signals and possible residual vari- 
ations by the variable comparison stars that could not be 
taken out. 

3.2.1 Comparison with literature data 

Heynderickx (1992) reported the detection of four frequen- 
cies in his light curves of V433 Car. Two of them agree with 
the two strongest signals in our Table 3, and the other two 



are 1 c/d aliases of our values. We prefer the frequencies 
derived our analysis because our data are more numerous, 
the results in the different filters were consistent and as our 
two-site measurements have a much better spectral window 
function. 

The order in which we detected the variation frequencies 
of V433 Car is different from the one obtained by Heynder- 
ickx (1992). Although this author does not list the ampli- 
tudes of the signals he determined, inspection of his pub- 
lished PDM periodograms suggests that the pulsations suf- 
fered amplitude variability in the 14 years that have elapsed 
between his and our observations. 



3.3 IL Vel 

The spectral window and amplitude spectra of our mea- 
surements of IL Vel are shown in Fig. 5. Two independent 
frequencies dominate the periodogram. After prewhitening 
those, a third signal stands out, and including that one in 
the multifrequency solution allows the detection of two com- 
binations frequencies. Combining the residuals in all filters, 
scaling the U data by a factor of 0.76 and the V data by 
1.05, results in the uninformative periodogram in the lowest 
panel of Fig. 5. The frequency solution for IL Vel is given in 
Table 4. 

This analysis may appear simple, but it is not. Firstly, 
the third independent signal is located halfway between the 
first two that are of much higher amplitude. Although it 
is well resolved from the two strong modes, some doubts 
remain as to its reality; it could be an artifact caused by 
amplitude/frequency variability of the two main signals. 

We have therefore divided our data set into two halves 
(before and after JD 2452360) and analysed them separately. 
We adopted the frequencies of the two strongest modes 
and their detected combination signals as definite and fit- 
ted them to each half, but left their amplitudes and phases 
as free parameters. We then computed amplitude spectra of 
the residuals for each half in each filter. The signal fs was 
detected in each of the 8 subsets of data; we conclude that 
it is real. 

The second cause of concern is the rms residual of our 
five-frequency fit to the IL Vel light curves. After combi- 
nation of the data in all four filters, it is an enormous 7.8 
mmag per point, twice as large as for V433 Car and 2.6 
times larger than for KZ Mus. The noise level in the lowest 
panel of Fig. 5 is also correspondingly higher. 
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Table 3. Multifrequency solution for our V433 Car data. Frequencies are mean values from results of all four filters weighted by the S/N 
ratio of the signal. Errors on the amplitudes were calculated with the formulae of Montgomery & O'Donoghue (1999). The S/N ratio 
quoted is the average of the values in the individual niters. 



ID 


Frequency 




Amplitude 




mean S/N 




(c/d) 


U (mmag) 


B (mmag) 


V (mmag) 


v (mmag) 




h 


9.12909 ± 0.00014 


9.9 ± 0.5 


8.2 ± 0.3 


8.2 ± 0.3 


8.8 ± 0.5 


17.9 


h 


8.31649 ± 0.00009 


10.0 ± 0.5 


8.2 ± 0.3 


7.3 ± 0.3 


8.0 ± 0.5 


16.7 


h 


9.66707 ± 0.00024 


4.3 ± 0.5 


3.8 ± 0.3 


4.0 ± 0.3 


4.2 ± 0.5 


8.5 


h 


7.7214 ± 0.0012 


3.4 ± 0.5 


2.8 ± 0.3 


2.3 ± 0.3 


2.6 ± 0.5 


5.4 


h - h 


0.53798 


1.8 ± 0.5 


1.2 ± 0.3 


1.5 ± 0.3 


1.3 ± 0.5 


3.7 



Table 4. Multifrequency solution for our IL Vel data. Frequencies are mean values from results of all four filters weighted by the S/N 
ratio of the signal. Errors on the amplitudes were calculated with the formulae of Montgomery & O'Donoghue (1999). The S/N ratio 
quoted is the average of the values in the individual niters. 





ID 


Frequency 
(c/d) 


U (mmag) 


Amplitude 
B (mmag) V (mmag) 


v (mmag) 


mean S/N 




h 


5.45976 ± 0.00012 


59.9 ± 0.8 


45.5 ± 0.6 


43.2 ± 0.7 


45.6 ± 0.7 


43.2 




h 


5.36293 ± 0.00011 


53.0 ± 0.8 


40.7 ± 0.6 


38.3 ± 0.7 


40.6 ± 0.7 


38.4 




h 


5.41340 ± 0.0010 


8.7 ± 0.8 


7.1 ± 0.6 


6.5 ± 0.7 


7.3 ± 0.7 


6.6 


h 


+ h 


10.82269 


5.2 ± 0.8 


4.0 ± 0.6 


4.1 ± 0.7 


3.7 ± 0.7 


4.2 


h 


-h 


0.09683 


5.3 ± 0.8 


3.3 ± 0.6 


3.4 ± 0.7 


4.4 ± 0.7 


3.7 



We started our quest for the reason for the high resid- 
uals by investigating the constancy of the comparison star. 
We calculated differential light curves to all the other stars 
that we had simultaneous data for, but found no evidence 
for variability of HD 79670, at least at a level that could 
produce the high residuals observed. They must therefore 
be (mostly) due to IL Vel. 

When we checked the reality of /3, we noted that its 
location in the amplitude spectra of the two subsets of data 
was somewhat shifted and that there was more structure 
surrounding it than just the normal spectral window pat- 
tern plus noise. A plot of the residual light curves after 
subtraction of our five- frequency solution shows conspicuous 
variability on time scales similar to those of the known pul- 
sational signals. We are therefore inclined to think that fur- 
ther, presently unresolved, pulsation frequencies are present 
in the light curves of IL Vel and that a data set with an even 
longer time base than ours is required to detect them. 



3.3.1 Re-analysis of published data 

The results of our frequency analysis differ from that by 
Heynderickx & Haug (1994). Our first frequency is signifi- 
cantly different from theirs, but we agree on the second one. 
Heynderickx & Haug (1994) did not find our third frequency, 
whereas we did not detect their third and fourth signals. 

To clarify the situation, we have retrieved the data by 
Heynderickx & Haug (1994) and re-analysed them. We show 
the periodogram analysis of their Johnson B data in Fig. 6, 
where we have assumed our frequencies to be correct, and 
used them as the initial values for light curve-fitting and 
prewhitening. 

Our two frequencies /i and fi are sufficient to ex- 
plain the dominant variations in these measurements. Af- 



ter they are prewhitened, the residual amplitude spectrum 
(lowest panel of Fig. 6) shows no peak exceeding 5 mmag; 
the two additional frequencies claimed by Heynderickx & 
Haug (1994) had B amplitudes of 20 and 14 mmag, respec- 
tively. We conclude that these were artifacts originating from 
choosing a wrong alias frequency. 

The absence of our signal fa in their data may be due 
to a similar reason. Two out of their three yearly data sets 
spanned two weeks, and one consisted two single-week runs 
spread by 25 days. The beating phenomenon of fi with ei- 
ther /i or f'2 is therefore never properly resolved. Due to 
this particular unfortunate time distribution fa can be sup- 
pressed below detectability, as we have checked with numer- 
ical simulations. It is therefore well possible that our fa was 
present in their data but could not be detected. We note 
that the two combination frequencies detected in our data 
may also be present in the measurements by Heynderickx & 
Haug (1994), but if so, their amplitude was lower. 

Finally, we would like to point out that the lowest panel 
of Fig. 6 also shows a residual mound of power, even more 
pronounced than in the residual periodogram of our data. 
The combined, scaled and average UBV residuals of the data 
by Heynderickx & Haug (1994) have an rms scatter of 9.8 
mmag, also much higher than their measurement accuracy. 
As these authors used a comparison star different from ours, 
we strengthen the conclusion that there are unresolved fre- 
quencies in the light curves of IL Vel. 

3.4 Variable comparison stars 

3.4.1 HD 90434 

This was the original local comparison star for V433 Car. 
When we examined amplitude spectra of the differential 
data reduced during the course of the observing campaign, 
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Figure 4. Spectral window and amplitude spectra of our B filter 
data of V433 Car. Some prcwhitening steps of the detected fre- 
quencies are shown in consecutive panels. The lowest panel shows 
the amplitude spectrum of the combined UvBV residuals after all 
previously significant signals have been removed from the mea- 
surements. One combination frequency is detected that way. 



we noted some additional low-frequency variability besides 
the pulsations of V433 Car. The differential light curves of 
HD 90434 relative to the comparison star for IL Vel, allowed 
us to identify HD 90434 unambiguously as a new variable. 
Figure 7 contains the spectral window and amplitude spec- 
tra of these data. 

We detect one variation frequency for this star. 
Prewhitening it from the data, combining the residual light 
curves as done for the (3 Cephei stars, and computing the 
amplitude spectrum of this data set results in some evi- 
dence for additional variability (lowest panel of Fig. 7) , but 
produces no more detections. The increase in noise level at 
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Figure 5. Spectral window and amplitude spectra of our B filter 
data of IL Vel. Some prewhitening steps of the detected frequen- 
cies are shown in consecutive panels. The lowest panel shows the 
amplitude spectrum of the combined UvBV residuals after all 
significant signals have been removed from the measurements. 



low frequencies may however also (partly) be due to residual 
effects of extinction and sky transparency as the two stars 
are separated by about 10° in the sky. The frequency solu- 
tion we derived for HD 90434 is given in Table 5. We note 
that the amplitude of the variability increases toward the 
blue, as confirmed by an examination for colour variations. 
No significant phase shifts between the light curves in the 
different filters were found. 



34.2 HD 89768 

After discovery of its variability, HD 90434 was replaced as a 
comparison star for V433 Car by HD 89768 whose spectral 
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Frequency (c/d) 

Figure 6. Amplitude spectra of Heynderickx & Haug's (1994) B 
filter data of IL Vel; consecutive panels show some prewhitening 
steps. The positions of our frequencies /i and are indicated 
with arrows and adopted as definite for the prewhitening. 




Frequency (c/d) 

Figure 7. Upper two panels: spectral window and amplitude 
spectrum of our B filter data of HD 90434 relative to HD 79670. 
Lower panel: prewhitened amplitude spectrum of the combined 
residual light curves out to the Nyquist frequency. 



Table 5 . Results of the frequency analysis for the two comparison 
stars. Frequencies are mean values from results of all four filters 
weighted by the S/N ratio of the signal. Errors on the amplitudes 
were calculated with the formulae of Montgomery & O'Donoghue 
(1999). The S/N ratio quoted is the average of the values in the 
individual filters. 



Star 


HD 90434 


HD 89768 


Frequency (c/d) 


1.5229 ± 0.0012 


0.6875 ± 0.0004 


U Amplitude (mmag) 


15.9 ± 1.2 


13.6 ± 1.2 


B Amplitude (mmag) 


11.5 ± 0.8 


14.8 ± 0.6 


V Amplitude (mmag) 


9.2 ± 0.7 


12.7 ± 0.5 


v Amplitude (mmag) 


11.4 ± 1.0 


15.5 ± 0.9 


Mean S/N 


5.9 


9.0 



type of Al IV/V places it outside any known pulsational 
instability region, but we were no luckier with this star. In 
fact, we only noticed that this object was variable as well 
when it was too late to replace it with yet another star. 

We again constructed differential light curves with re- 
spect to HD 79670 and show its Fourier analysis in the same 
fashion as the one in the previous section in Fig. 8. Again, a 
single frequency suffices to explain the observed light vari- 
ations, and some residual low-frequency variability seems 
present that we attribute to atmospheric effects. Table 5 
contains our results for HD 89768 as well. 

From the amplitudes in the different filters listed in Ta- 
ble 5 it appears that this star also shows some colour vari- 
ability. We have therefore analysed the colour light curves as 
well, but we could not confirm this suspicion. It seems that 
the different amplitudes in the four filters are generated by 
noise effects and our errors are underestimates for variations 
at these low frequencies. 



4 BASIC STELLAR PARAMETERS 
4.1 The (3 Cephei stars 

The astrophysical analysis of our results is facilitated 
by firstly restricting the physical parameter space oc- 
cupied by these stars. Consequently, we have estimated 
their masses and temperatures with calibrations of mul- 
ticolour photometry, as meaningful trigonometric paral- 
laxes are not available. Fortunately, all three stars have 
already been measured in the Str0mgren, Geneva and 
Johnson systems. We have retrieved the corresponding 
standard values from the Lausanne-Geneva data base 
(http : //obswww.unige . ch/gcpd/gcpd.html ). 

We have then applied different calibrations for effective 
temperature, surface gravity and absolute magnitude to the 
three stars. It seems important to use more than one cali- 
bration as their results often deviate. We can then hope to 
get a good average result with a realistic determination of 
its uncertainty. Absolute visual magnitudes were calculated 
from the Str0mgren (3 index using Crawford's (1978) results 
that also allow a determination of the interstellar reddening. 
Wc find A v = 0.90 for IL Vel, A v = 0.38 for V433 Car and 
A v = 1.24 for KZ Mus. 

We used the formulae by Napiwotzki, Schonberner & 
Wenske (1993) to calculate stellar temperatures from the 
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Figure 8. Upper two panels: spectral window and amplitude 
spectrum of our B filter data of HD 89768 relative to HD 79670. 
Lower panel: prewhitened amplitude spectrum of the combined 
residual light curves out to the Nyquist frequency. 



Table 6. The adopted effective temperatures and luminosities for 
our target stars from standard photometry. 



Parameter 


IL Vel 


V433 Car 


KZ Mus 


T ef f [kK] 


23.6 ± 0.6 


26.6 ± 0.7 


26.0 ± 0.7 


log L/Lq 


4.19 ± 0.22 


4.20 ± 0.20 


4.22 ± 0.20 



Str0mgren [u — b] and (b — y)o indices as well as from John- 
son (B — V)o. The tables by Flower (1996) allow another 
temperature estimate from (B — V)o and also provide bolo- 
metric corrections. The latter were also determined from 
Drilling & Landolt's (2000) tables. Estimates of the surface 
gravities of the three stars were obtained with the calibra- 
tion by Smalley & Dworctsky (1995) from the Str0mgren f) 
index, including a correction proposed by Dziembowski & 
Jerzykiewicz (1999). Finally, we applied the work by Kiinzli 
et al. (1997) for the Geneva system that results in effective 
temperatures and surface gravities. Average values for the 
resulting effective temperatures and luminosities as well as 
error estimates are given in Table 6. 



4.2 The variable comparison stars 

The HIPPARCOS parallax of HD 90434 (ESA 1997) has a 
relative error that is too large to make it useful. The spectral 
classification and the Str0mgren colours of the star imply it 
is a late B star. Johnson and Geneva colours of this object 
are also available at the Lausanne-Geneva data base. Con- 



sequently, we applied the same calibrations that we have al- 
ready used for the /3 Cephei stars to determine its effective 
temperature and luminosity. We derive T e g = 12.2 ± 0.4 kK 
and log L/Lq = 1.86 ± 0.10. 

No standard photometry is available for HD 89768. 
However, from the mean magnitude differences of V433 
Car from both of its comparison stars, we can estimate 
(B — V) — 0.09 and (U - B) = 0.04 for HD 89768. With 
the HIPPARCOS parallax of 2.09 ± 0.69 mas and the galac- 
tic reddening law by Chen et al. (1998), we can then es- 
timate (B - V) = 0.01 ± 0.03 and M v = -0.7 ± 0.9. 
From the calibration by Napiwotzki et al. (1993) we then 
obtain T g — 9250 ± 400 K, and the luminosity results in log 
L/L Q = 2.3 ±0.4. 

4.2.1 The nature of HD 90434 and HD 89768 

The effective temperature and luminosity of HD 90434 as 
derived above put it into the instability region of the slowly 
pulsating B (SPB) stars (Dziembowski, Moskalik & Pamyat- 
nykh 1993). The colour variability we found combined with 
the variation being in phase in the different filters supports 
this hypothesis and argues against an interpretation in terms 
of rotational modulation of a chemically peculiar star; the 
Str0mgren mi index for HD 90434 does not suggest it is an 
Ap star either. 

On the other hand, HD 89768 is not located in any 
known pulsational instability region, and we could not detect 
colour variability in its light curves. We suspect that the 
variability may be binary-induced, but spectroscopic data 
would be required to confirm our hypothesis. 



5 MODE IDENTIFICATION 

The ranges in effective temperature and luminosity of the 
(5 Cephei stars listed in Table 6 are an essential ingredi- 
ent for mode identification. We have followed the method 
proposed by Balona & Evers (1999, see their paper for a de- 
tailed description). It uses theoretically calculated nonadia- 
batic parameters to determine the amplitude ratios between 
different wavebands that are the discriminator between the 
different modes of pulsation. 

We computed stellar evolutionary models by means of 
the Warsaw-New Jersey evolution and pulsation code (de- 
scribed, for instance, by Pamyatnykh et al. 1998) for solar 
chemical composition. The models spanned a mass range 
between 8 — 16Mq in steps of IMq. Guided by the pro- 
jected rotational velocities of the stars listed in Table 1, 
we chose a rotational velocity of 100 km/s on the ZAMS 
for models representing IL Vel and KZ Mus, but 260 km/s 
for V433 Car. By comparing the model evolutionary tracks 
with the stellar parameters in Table 6, we inferred that IL 
Vel is a star with a mass between 10.5 — 13Mq in the sec- 
ond half of its main-sequence evolution, whereas KZ Mus 
and V433 Car are less evolved but somewhat more massive 
(11.5M Q < M, < 14M Q ). 

We then proceeded by calculating theoretical UvBV 
amplitudes for all pulsationally unstable modes along model 
sequences that spanned the physical parameter space of our 
targets, whereby we extended the error bars in Table 6 by 
up to a factor of 2 to be conservative. We did not consider 
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the phase shifts between the different photometric bands 
because their measured relative error is much higher than 
that for amplitude ratios. We have not found phase shifts 
significant at the 3a level and they have therefore no dis- 
criminative power. 

The resulting theoretical amplitudes in the four pass- 
bands were normalised to that in the U filter, and the av- 
erage amplitude ratio and its rms error were calculated. In 
this way we can also get an idea of the uncertainties in- 
volved in the calculation of theoretical results that need to 
be considered when being matched with the observations. 

The comparison between the theoretical and observed 
colour amplitude ratios of IL Vel is shown in Fig. 9. We also 
include the UBV colour amplitude ratios from the data by 
Heynderickx & Haug (1984) with the correct frequencies /i 
and / 2 . 

The identification for the two strongest modes of IL 
Vel is unambiguous: both have a spherical degree 1=1. 
The third pulsation mode we found that has a frequency 
intermediate between the other two is also nonradial and 
could be I = 1,2 or 3. Our mode identification is therefore 
consistent with that by Heynderickx et al. (1994), although 
they have used an incorrect alias and spurious frequencies. 

Observed and theoretical amplitude ratios for the four 
modes of V433 Car are plotted in Fig. 10. The identifications 
we derive from that figure are not as clear as the ones for IL 
Vel. Mode /i is most likely £ = 2, and the amplitude ratios 
for / 2 agree best with I — 1. All we can say about mode /3 
is that it is nonradial (t — 2 or 4?), and mode fi may be 
I = 0, 1, 2 or 3. 

Heynderickx et al. (1994) arrived at similar conclusions 
concerning the identification of the modes of V433 Car, and 
they also reported inconsistencies of their identifications in 
different photometric systems. We must therefore consider 
the attempts to derive mode identifications for V433 Car as 
mostly inconclusive. 

Turning now to KZ Mus, we show the comparison of 
its observed amplitude ratios with their theoretical predic- 
tions in Fig. 11. The identifications for the highest-amplitude 
modes of this star are very clear. The strongest mode of KZ 
Mus is quadrupole (£ = 2), whereas mode / 2 is radial. This 
is in excellent agreement with the identification by Aerts 
(2000). 

As we have a much larger data set available than Aerts 
(2000) did, we can also provide an identification for the third 
mode of KZ Mus, which is a dipole {£= 1). Only the identi- 
fication for the weak fourth mode of this star does not point 
to an unambiguous spherical degree, but we can safely that 
it is nonradial with I < 3; a radial identification is ruled 
out because this mode's frequency is too close to that of the 
radial mode / 2 . 

Finally, we also look for a possible mode identification 
for our SPB candidate HD 90434. For this star, we adopted 
3Mq models near the ZAMS guided by the results from 
Sect. 4.2 to determine the theoretical UvBV amplitude ra- 
tios. We then proceeded in the same manner as for the 13 
Cephei stars, and present the comparison between theoreti- 
cal and observed amplitude ratios in Fig. 12. 

Modes of spherical degrees I — 1 and 4 give the best 
agreement, but 1 = 2 also seems a viable identification. The 
long period of the star and the normalisation of the ampli- 
tudes in the U band (where the absolute error on the ampli- 
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Figure 9. Observed and theoretical UvBV amplitude ratios 
(lines) for < I < 4 for IL Vel. Amplitudes are normalised at 
U. The filled circles with error bars are from our measurements, 
the open circles (somewhat shifted in wavelength for reasons of 
display) with error bars are for the data by Heynderickx & Haug 
(1984). The full lines are for radial modes, the dashed lines for 
dipole modes, the dashed-dotted lines for quadrupole modes, the 
dotted lines for modes of i = 3 and the dashed-dot-dot-dotted 
lines are for £ = 4. The small error bars, also slightly shifted in 
wavelength, denote the uncertainty in the theoretical amplitude 
ratios. The upper panel is for mode /1, the middle one for / 2 , and 
the lower one for fy. 



tude is largest) are suspected to result in larger systematic 
contributions to the observational errors than the formal 
values we adopted would include. Taking the effects of geo- 
metrical cancellation (Dziembowski 1977) into account, the 
most likely reason for the variability of HD 90434 is pulsa- 
tion in a dipole mode. 



6 ASTEROSEISMOLOGY 
6.1 IL Vel 

We will now take advantage of the information gathered 
in the previous sections to understand the pulsational be- 
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Figure 10. Observed (filled circles with error bars) and theo- 
retical UvBV amplitude ratios (lines) for < I < 4 for V433 
Car. Amplitudes are normalised at U. The full lines are for radial 
modes, the dashed lines for dipole modes, the dashed-dotted lines 
for quadrupole modes, the dotted lines for modes of £ = 3 and 
the dashed-dot-dot-dotted lines are for £ = 4. The small error 
bars somewhat shifted in wavelength denote the uncertainty in 
the theoretical amplitude ratios. The upper panel is for mode fi , 
the second for /2 , the third for , and the lower one for fn . 
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Figure 11. Observed (filled circles with error bars) and theoreti- 
cal UvBV amplitude ratios (lines) for < £ < 4 for KZ Mus. Am- 
plitudes are normalised at U. The full lines are for radial modes, 
the dashed lines for dipole modes, the dashed-dotted lines for 
quadrupole modes, the dotted lines for modes of £ = 3 and the 
dashed-dot-dot-dotted lines are for £ = 4. The small error bars 
somewhat shifted in wavelength denote the uncertainty in the 
theoretical amplitude ratios. The upper panel is for mode /i, the 
second for /2, the third for /3, and the lower one for fi. 
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Figure 12. Observed (filled circles with error bars) and theoret- 
ical UvBV amplitude ratios (lines) for 1 < t < 4 for HD 90434. 
Amplitudes arc normalised at U. The dashed lines are for dipole 
modes, the dashed-dotted lines for quadrupole modes, the dotted 
lines for modes of t = 3 and the dashed-dot-dot-dotted lines are 
for t = 4. The small error bars somewhat shifted in wavelength 
denote the uncertainty in the theoretical amplitude ratios. 



haviour of our three target stars by means of comparison 
with model calculations. We have again used the Warsaw- 
New Jersey stellar evolution and pulsation code for this part 
of the analysis. We note that we do not attempt an exhaus- 
tive seismic analysis of the stars because too few identified 
pulsation modes are available for such purposes. Instead, we 
want to determine what we can learn about the stars by 
means of standard models and procedures and we want to 
assess the potential for asteroseismology of /3 Cephei stars. 

One tool that has been established for constraining the 
range of physical parameters of the seismological model of a 
star is stability analysis (see Pamyatnykh 2003 for a review). 
As a model of a pulsating star evolves, its radius changes, 
and the range of pulsationally unstable radial overtones of 
the eigenfrequencies varies as well. In conclusion, the fre- 
quency range of the unstable eigenmodes depends on the 
location of the model in the instability strip. This frequency 
range can then be matched with the one actually observed 
in the star. 

Applying this method to IL Vel, we can take advantage 
of our mode identification. As all the modes for which we 
have a certain identification are I = 1, we can restrict the 
comparison to this single spherical degree only. Figure 13 
shows evolutionary tracks for models between 8 and 16 Mq 
(we assumed solar chemical composition and a rotational 
velocity of 100 km/s on the ZAMS) in the theoretical HR 
diagram together with the parameters of the star from Table 
6; models with a matching range of unstable frequencies are 
indicated. 

Inspection of Fig. 13 shows very good agreement be- 
tween the star's position in the HRD and the range of mod- 
els that have the same frequency domain unstable. We note 
that the sequences of matching models sometimes have gaps 
as they evolve along the main sequence. This is due to a 
change in the type of modes that produce the match and 
their frequency evolution which is affected by avoided cross- 
ings. We caution that the strength of pulsational driving in 
(5 Cephei stars is strongly dependent on metal abundance 
(Moskalik & Dziembowski 1992, Pamyatnykh 1999). There- 
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Figure 13. The position of IL Vel in the theoretical HR dia- 
gram. Some stellar model evolutionary tracks labeled with their 
masses (full lines) and the theoretical borders of the (i Cephei star 
instability strip (Pamyatnykh 1999, dashed lines) are included 
for comparison. The sections of the evolutionary tracks that are 
marked with filled circles represent models in which the observed 
frequency range of the star is unstable. 



fore the results in Fig. 13 are only an indication of the range 
of possible models for IL Vel. 

The two £ = 1 modes of IL Vel have similar frequencies; 
they are spaced by only 0.0986 c/d. The projected rotational 
velocity of the star (wsini = 65 km/s, cf. Table 1) implies a 
rotational period P ro t < 5.8 d for the star assuming a mass 
of 12 Mq as indicated by Fig. 13. At first glance, this sug- 
gests that those two modes cannot be caused by rotational 
splitting of a single mode. Three explanations for this close 
splitting may be possible. First, the star's pulsations are just 
undergoing an avoided crossing phenomenon (e.g. see Aizen- 
man, Smeyers & Weigert 1976) resulting in a close proximity 
of a p- and a g-mode. Second, the two I = 1 modes are the 
central and prograde components of a multiplet, and the 
second-order rotational splitting is large enough to produce 
such a small frequency difference. Third, we may observe a 
gravity mode; the first-order rotational splitting of such a 
mode is approximately v £s (1 — (£(£ + l)) _1 )P rot 1 (Winget 
et al. 1991). 

We examine these hypotheses in Fig. 14, where we show 
the evolution of the rotational m-mode splitting and the 
frequencies of the I = 1 modes of a 12 Mq model with 
frot.zAMS = 70 km/s along the main sequence. The first 
three pressure (p) and the first two gravity (g) modes are 
shown. From the upper panel of this figure it becomes clear 
that the observed frequency splitting of IL Vel can be easily 
reproduced by the model in the corresponding temperature 
range. In fact, this would be possible with models rotating 
as fast as 110 km/s. On the other hand, the rotational fre- 
quency splitting of any of these modes never becomes small 
enough that the observed frequency difference can be due to 
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Figure 14. Upper panel: the size of rotational splitting for t = 
1 modes of a 12 Mq model with « ro t,zAMS = 70km/s along 
its evolution along the main sequence. The first three p-modes 
and the first g-mode are indicated. The horizontal dashed lines 
denote the frequency difference between the two strongest modes 
of IL Vel, and the vertical solid line with parallel dotted lines 
is the effective temperature of the star and its error estimate 
from Table 6. Lower panel: stability and frequency evolution of 
t = 1, m = mode frequencies in the same model. Modes that 
are pulsationally unstable are denoted with filled circles, the open 
circles are for stable modes. The two dashed lines represent the 
observed 1 = 1 frequencies, and the vertical lines again labels the 
effective temperature of the star and its error estimate. Note the 
avoided crossing of the pi and the gi mode starting at about log 
T efI = 4.39. 



£ = 1 modes of the same type and radial overtone but with 
m = 1 and m = — 1. This also implies that the third mode 
of IL Vel cannot be another member of this multiplet. 

To determine the modes that are potentially responsible 
for the pulsations of IL Vel, we have traced the frequencies of 
the individual I = 1 modes of the same model as used before 
along the main sequence in the lower panel of Fig. 14. The 
mode originating as pi on the ZAMS has a frequency very 
similar to the two observed for IL Vel over the whole temper- 
ature range estimated for the star, and is unstable as well. 
This mode has already undergone an avoided crossing with 
the gi mode and should therefore also have gravity mode 
characteristics, i.e. it is a mixed mode. In higher-mass mod- 
els, dipole modes at this frequency become stable against 
pulsations, and in models below IIMq the gi mode re- 
produces the observed frequency in the temperature range 
appropriate for IL Vel. 

Another observation from Fig. 14 is that the pi and g\ 
modes never come closer together than about 0.5 c/d. It is 
therefore unlikely that the two modes observed in IL Vel 
originate from different modes that just happen to undergo 
an avoided crossing. If we consider the hypothesis that the 
proximity of the two observed I = 1 modes is caused by 
the effects of second-order rotational splitting, we find that 
frot ~ 160 km/s, and hence i < 25°, which is also unlikely. 
In addition, the high amplitude of the prograde mode and 
the absence of the retrograde component of the hypothesised 
triplet would also be difficult to explain. 

We conclude that the two dipole pulsation modes of IL 
Vel are most likely due to rotationally split components of 
the mode originating as p\ on the ZAMS. One of these modes 
is m = 0, the other one \m\ = 1, and 65 km/s < v rot < 110 
km/s. 

6.2 V433 Car 

We have performed a similar analysis for V433 Car. Regret- 
tably, we could not obtain unambiguous £ identifications for 
any of its modes. On the other hand, the excited frequency 
range of the star is larger, which may however be related 
to it being a fast rotator. In any case, we show the match 
of unstable model frequencies to that of the observed ones 
in Fig. 15. The model evolutionary tracks are again for so- 
lar chemical composition, but for a rotational velocity of 
260 km/s on the ZAMS, and we considered all modes with 
< £ < 4. 

Again, the agreement between theory and observation 
is fairly good, even when keeping in mind that our analysis 
does not take the frequency spread due to fast rotation of 
the star. However, as the observational results only allow 
us a rough qualitative comparison between observations and 
theory for V433 Car, this is not regarded as being important. 
V433 Car appears to be an object of about 13 Mq that has 
just entered the (3 Cephei star instability strip. 

6.3 KZ Mus 

The asteroseismological prospects of KZ Mus appear much 
better. In particular, we have unambiguously identified one 
of its modes as radial, and two more modes also have unique 
£ assignments. In the same fashion as in the previous sec- 
tions, we show model evolutionary tracks (for ?) ro t,zAMS = 
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Figure 15. The position of V433 Car in the theoretical HR dia- 
gram. Some stellar model evolutionary tracks labeled with their 
masses (full lines) and the theoretical borders of the {3 Cephei star 
instability strip (Pamyatnykh 1999, dashed lines) are included 
for comparison. The sections of the evolutionary tracks that are 
marked with filled circles represent models in which the observed 
frequency range of the star is unstable. For some models, we ob- 
tain unstable pulsation modes somewhat outside the blue edge of 
the theoretical instability strip, which is due to I = 4 modes; the 
borders indicated are only for I < 2. 

100 km/s) in the theoretical HR diagram, the star's position 
in it, and models with unstable modes in the observed fre- 
quency range in Fig. 16. In addition, we have indicated the 
locations of models that have a radial mode at 5.9506 c/d 
as observed for KZ Mus. 

As for the other two stars, we find no disagreement 
between models with unstable modes in the observed fre- 
quency range and the observationally determined position 
of the star in the HRD. Figure 16 also allows us to reject 
the second overtone as the identification for the radial mode 
of KZ Mus, as it is inconsistent with the star's evolutionary 
state, and it never becomes pulsationally unstable. 

However, we cannot say at this point whether the radial 
mode is the fundamental or first overtone. Models with the 
fundamental radial mode at 5.9506 c/d are in better agree- 
ment with the stellar parameters from Table 6, but it may 
be stable in models with the mass inferred for the star. On 
the other hand, for this mode to be the first overtone, only 
models at the upper limits on mass and luminosity seem 
appropriate. 

We attempted to find further constraints on the pa- 
rameter space occupied by KZ Mus by using the nonradial 
modes that we identified unambiguously; any model repre- 
senting the star should have modes of the same £ in the 
correct frequency range. All models with the radial funda- 
mental or first overtone mode at 5.9506 c/d that we inves- 
tigated could be used to reproduce the observed nonradial 
modes given the unknown angular rotational velocity of the 
star. However, when imposing mode stability as an addi- 
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Figure 16. The position of KZ Mus in the theoretical HR dia- 
gram. Some stellar evolutionary tracks labeled with their masses 
(full lines) and the theoretical borders of the Cephei star insta- 
bility strip (Pamyatnykh 1999, dashed lines) arc included for com- 
parison. The sections of the evolutionary tracks that are marked 
with filled circles represent models in which the observed fre- 
quency range of the star is unstable. The lines going from the 
top left to the bottom right are locations of models that have a 
radial mode at 5.9506 c/d. For models on the dashed-dotted line 
this would be the fundamental radial mode, those connected by 
the dotted line show the first overtone at this frequency, and the 
ones connected by the dashed-dot-dot-dotted line have their sec- 
ond radial overtone at 5.9506 c/d. The thick parts of these lines 
connect the models which are pulsationally unstable. 

tional criterion, it is found that model sequences of 12 and 
13 Mq give best agreement with the observations. 

In this case, it is most likely that /2 is the radial fun- 
damental mode of KZ Mus. Then the £ — 1 mode would 
be associated with the mode originating as pi at the ZAMS 
and the observed 1 = 2 most likely corresponds to the mode 
originating as g\ . 



7 SUMMARY, DISCUSSION AND 
CONCLUSIONS 

In an attempt to investigate the asteroseismological poten- 
tial of the f3 Cephei stars, we have obtained 127 - 150 h of 
differential photoelectric photometry for each of the three 
pulsators IL Vel, V433 Car and KZ Mus. In line with the 
recent results by Stankov et al. (2002) and Cuypers et al. 
(2002), we have shown that low-amplitude pulsation modes 
are present in /3 Cephei stars, and that large observational 
efforts to determine their mode spectra in detail are justified. 
We therefore confirm the importance of (3 Cephei stars as 
asteroseismological targets for both ground and space based 
measurements. 

By using the amplitude ratios of the different pulsation 
modes in the UvBV filters of the Johnson and Str0mgren 
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systems, we performed mode identification. We obtained un- 
ambiguous results for all pulsation modes with photometric 
amplitudes larger than about 10 mmag in V, or with relative 
errors on their amplitude smaller than 3%. The modes we 
could identify were the two dominant ones of IL Vel and the 
three strongest of KZ Mus, but none of V433 Car. 

This may not only be due to its photometric amplitudes 
that are small compared to those of the other two (3 Cephei 
stars. We have shown that V433 Car is a fast rotator, and 
therefore mode coupling (see Daszynska-Daszkiewicz et al. 
2002 and references therein) may occur. Mode coupling af- 
fects the predicted colour amplitude ratios and phase shifts. 
Since we have not taken this into account in our analysis, it 
may be at least partly responsible for the ambiguous mode 
identifications for this star. 

On the other hand, our results for the more slowly ro- 
tating stars appear reliable, which represents a good motiva- 
tion for studying a larger number of multiperiodic /3 Cephei 
stars in the future. Our limit of the relative error of 3% 
of the photometric amplitudes for successful mode identi- 
fication can then serve as a guideline for planning future 
campaigns. 

We have performed model calculations for our three tar- 
gets. The outcome was also quite encouraging. Although we 
only detected up to four modes per star (thus our oberva- 
tional constraints on the models are not very strong), in no 
case disagreements between the observational results and 
model predictions were found. The observationally deter- 
mined positions of the stars in the HR diagram could be 
reproduced by models that have the same range in pulsa- 
tion frequencies unstable as are excited in the real stars. 

In addition, we were able to place constraints on the 
types of mode observed in IL Vel and KZ Mus. We have 
shown that the most likely explanation for the two dipole 
modes of IL Vel are rotationally split components (one nec- 
essarily being m = 0) of the mode originating as pi on the 
ZAMS. The third mode of IL Vel cannot be another multi- 
plet member, and the equatorial rotational velocity of the 
star must be smaller than 110 km/s. If the radial mode 
observed in KZ Mus is the fundamental, then the high- 
amplitude 1 = 2 mode is most likely the one originating 
as gi, and the I = 1 mode we identified is probably the one 
originating as pi . 

We cannot push the asteroseismological analysis of the 
three stars further at this point because too few modes have 
been detected and identified for detailed modelling. How- 
ever, we strongly suspect that more modes are present in 
IL Vel that just need better time resolution to be detected. 
KZ Mus may have more modes as well, and a unique mode 
identification for its fourth mode may place a constraint on 
the star's rotation frequency. Further work on both stars 
should also include time-resolved spectroscopic observations 
to determine the azimuthal order m of the high-amplitude 
nonradial modes; although both objects are V « 9.1, the 
spectra we obtained suggests this is possible. Determina- 
tions of the metal abundances of the stars would also be 
a valuable ingredient for model calculations. With all these 
observations in hand, precision asteroseismological studies 
of IL Vel and KZ Mus appear possible. 
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